IMSc/2002/ 12/43 
TIFR/TH/02-36 
December, 2002 



?7 C production and dimuon enhancement 
in heavy ion collisions 



Rahul Basu 1 

The Institute of Mathematical Sciences, 
C.I.T. Campus, Chennai 600 113, India. 

K. Sridhar 2 

Department of Theoretical Physics, 
Tata Institute of Fundamental Research, 
Homi Bhabha Road, Mumbai 400005, India. 



Abstract 

Dilepton production in heavy ion collisions, in the Intermediate Mass Region 
(IMR) has consistently shown an excess over theoretical estimates. An attempt 
to understand this discrepancy between the observed dilepton pairs and the 
theoretical estimate is made here through the production of the r/ c meson and 
estimates obtained by NRQCD calculations. We find that ?? c production offers 
a satisfactory quantitative picture for explaining the discrepancy. 
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Dilepton production plays a very important role in the study and understand- 
ing of heavy ion collisions. This is mainly because dileptons do not interact with 
the surrounding hadronic medium after being produced in a nucleus-nucleus col- 
lision. The typical dilepton invariant mass spectrum appears as a wide continuum 
interrupted by various resonance (p, u, 4>, J/tp, ip' etc) decay peaks. 

The Intermediate Mass Region (between the 4> and J/ip peaks, 1.5 GeV to 2.5 
GeV) is particularly interesting because it is believed to contain dileptons created 
in the thermalised QGP produced in nucleus-nucleus collisions (thermal dimuons). 
However it is precisely in this region that many different experiments in A-A col- 
lisions |1] have shown an excess of dimuon production. In all these data sets, the 
dilepton sources are either Drell-Yan pairs or decays of J/ip or DD. While the 
bulk of the data agrees with such a production picture, there is a significant dis- 
crepancy between the observed dilepton pairs and the theoretical estimate based 
on the above sources in this Intermediate Mass region (IMR) with p + pr invariant 
mass in the range 1.5 - 2.5 GeV. 

Many explanations have been offered in the literature for this excess viz. de- 
crease in p meson mass due to thermal effects in e + e~ data |2 J, D-rescattering 151 , 
enhanced DD production, in-flight 7r + 7r decaying to e + e~ |4J, fireball hydrody- 
namics |5J and so on. 

In the entire kinematic range (upto 5 GeV) other than the IMR region men- 
tioned above, where theory and experiment agree, the overall picture involving 
charm quarks |6| is that when there is sufficient energy exchange in a collision, 
protons have non-negligible charm content (cc pairs), and substantial high energy 
gluons which in turn can decay to cc pairs. These cc pairs occassionally form bound 
states such as J / ip by emitting a soft gluon to maintain color balance, or can further 
polarize uu or c/Jfrom the surrounding medium to form DD pairs. Although the 
present theoretical understanding cannot predict absolute numbers for these pro- 
cesses it is possible to check the consistency of this picture with various p-A and 
A-A data. By and large the data agrees with various quantitative checks. 

Other charm meson bound states can also be produced, such as rj c , ip' and x"s, 
though their relative abundance is constrained by size, the larger ones being less 
likely to be formed. 

However the r] c meson, which is a IS orbital state is expected to have the same 
size and has almost the same mass as the J /if). They only differ in their spin and 
hence in any collision where cc quarks are produced these can form rj c with about 
1/3 probability as J/ip. However the possibility of the production of r\ c mesons in 
nucleus-nucleus collisions was only considered very recently in wherein this 
was presented as a possible explanation of the discrepancy between theory and 
experiment in nucleus-nucleus collision - in particular in S-U and Pb-Pb collisions. 
However this paper presented a somewhat qualitatitive field theoretical picture 
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of the production of r] c mesons and used it to explain the discrepancy in the IMR 
region. Furthermore, this simple field theoretical picture was unable to account for 
the discrepancy seen in the IMR region for Pb-Pb collisions in the central region. 

In what follows we attempt to understand the discrepancy in the IMR region 
once again through the production of the r) c but by using a more quantitative de- 
scription for the production of the rj c meson. Our attempt may therefore, be seen 
as a way to build upon the ideas of Q using the technology of NRQCD to present 
a more quantitative picture. 

In NRQCD 13, the quarkonium production cross-section factorises into a 
perturb atively calculable short-distance (< I/Mq, Mq is the mass of the heavy 
quark) effect and a long-distance part which is given by non-perturbative matrix 
elements. The cross section for production of a quarkonium state H can be written 
as 



The coefficients F n correspond to the production of QQ in the angular momentum 
and colour state (singlet or octet) denoted by n and is calculated using perturba- 
tive QCD. The non-perturbative part, (0%) of mass dimension d n in NRQCD has a 
well-defined operator definition and is universal. These matrix elements can be ex- 
tracted from any one process and then be used to predict other processes where the 
same matrix elements appear. Though the summation involves an infinite num- 
ber of terms, the relative magnitude of the various terms is predicted by NRQCD 
and these matrix elements scale as powers of the relative velocity v. However, this 
does not necessarily imply that effects from higher orders in v will always be small 
in physical processes, because any observable, like the decay width or the cross 
section is given by a double expansion in the strong coupling constant « s (Mq) and 
the relative velocity v. 

The non-perturbative matrix elements in NRQCD, are not calculable and have 
to be obtained by fitting to available data. The matrix elements of the colour-singlet 
operators can be obtained from the quarkonium decay widths and the colour-octet 
matrix elements have been obtained by fitting NRQCD predictions to the CDF data 
11101 (111. The remarkable thing is that the non-perturbative parameters appearing 
in the r) c production cross section can be determined from the matrix elements de- 
termined from J / ip production at the Tevatron: this happens because of the heavy- 
quark symmetry of the NRQCD Lagrangian. This has been exploited earlier in the 
context of h c and r\ c production at Tevatron lTT2l[T3l. 

A Fock space expansion of the physical r\ c , which is a 1 S (J PC = 0~ + ) state 
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yields 

\ Vc ) = 0(1) QQ[ l S [ o ] }) + 0{v 2 ) QQ[ l Pl 8] ] 9) + 0(v 4 ) QQ[ 3 S [ ?] g) + ■ ■ ■ . (2) 



The colour-singlet 1 S state contributes at (9(1) but the colour-octet l P\ and 3 S'i 
channels effectively contribute at the same order because the P-state production is 
itself down by a factor of 0(v 2 ). The colour-octet states become a physical r/ c by 
the 1 P 1 ' 8 ' state emitting a gluon in an El transition, and by the 3 sj 8 ' state emitting a 
gluon in an Ml transition. The contributing subprocess cross sections are 



qq -> QQ[ 2S+1 l 



9 9 - QQ[ Lj], 

where the is in the 1 S [ o\ 1 S l o ] and 3 sf ] states. The ^f 1 state does not contribute 
at this leading order in a s because of Yang's theorem. 

We have computed the contributions to the cross-section for rj c production from 
the x sf, x sf and 3 sf ] states. 

Heavy quark spin-symmetry is made use of in obtaining (O^) 's from the exper- 
imentally available (O^J^)'s. Using this symmetry we get the following relations 
among (0*)'s: 

(OlOff^ollO) = (0|O^[ 3 5i]|0> (l + 0(v 2 )), 

(0| O^fS,} |0) = (0| O^So] |0) (1 + 0(v 2 )). (3) 

For the singlet matrix elements we have (0| 0^[ 3 Si] |0) = 1.2 GeV 3 The CDF J/ip 
data only constrains a combination of octet matrix elements given by A\ + A 2 = 

<°l°^l°> + = (2 . 2 ±0 . 5) x 10 -2 GeV s m CDF J/v , data donot 

allow for a separate determination of the values of A\ and A 2 because the shapes 
of these two contributions to the J/ ip pr distribution are almost identical. For our 
numerical predictions we assume that the non-perturbative matrix element of in- 
terest to us ((0| Oq c [ 3 Si] |0)) lies in the range determined by this sum. 

Using these values for the non-perturbative matrix elements we can compute 
the cross-section for rj c production. This has to be convoluted with the branching 
ratio for the r\ c — > 77* — > 7/i+zi - . This quantity can be estimated only approx- 
imately and to do this we have used the value of the r/ c — > 77 branching ratio 
modulated by a (1 — M 2 */M 2 J in the denominator. With these inputs, we have 
computed the /i + /i~7 yield in Pb-Pb and S-U collisions. The nuclear parton dis- 
tributions have been taken from the parameterisation of EKS |14|. The cuts and 
acceptances for the muons and photons have been taken from the experimental 
papers (see, for example, I15ID . 
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Figure 1: The graphs show the S-U (a) and Pb-Pb (b) results for central collisions. 
Note that in the IMR region the rj c contribution appears to saturate the deficit. 



Our results for dN/ dM are compared to the experimental curves for central col- 
lisions in in Fig. 1. In keeping with the convention used by the experimentalists, 
we have indicated the r\ c production contribution separately (long-dashed line) as 
well as the summed contribution of all the individual contributions in the IMR re- 
gion (dashed-dotted line). It appears from both the graphs that the r\ c contribution 
saturates the discrepancy seen in the IMR region. However this agreement is based 
on certain assumptions. First of all, the branching ratio of r\ c — ► 77* — > 7/i + /i~ 
used, as already stated in the previous paragraph, is only approximately estimated 
through very general field theoretical arguments. Secondly, as with all leading or- 
der calculations, the result is dependent (though mildly) on the scale of a s used in 
the calculation. Finally since, as we have explained, the CDF J /if) data only con- 
strains the combination of octet matrix elements given by A\ and A 2 we have, for 
convenience, assumed that the non-perturbative matrix element < 0|d?g / '^[ 1 <S'o]|0 > 
(i.e A 2 ) saturates the sum and hence gives the value of the < 0|Og e [ 3 5'i]|0 > matrix 
element. However this turns out to be not such a serious approximation and we 
find that varying this value between A\ and A 2 produces very little change in the 
overall shape and magnitude of the i] c contribution. 

The general arguments for the rj c contribution used in |7[ were unable to explain 
the discrepancy in central collision in Pb-Pb collision and it was speculated there 
that the difference between theory and experiment even after including r\ c produc- 
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tion, could perhaps be accounted for by glueball production. It seems however, in 
the present, more quantitative, analysis, that this possibility is ruled out since the 
dileptons from i] c seem to cover the deficit. In any case, even if such a scenario of 
glueball production is envisaged, the contribution would clearly be very small. 

Finally we have made no comments about peripheral collisions where also this 
discrepancy in the IMR region is seen. This is because a quantitative study of 
peripheral collision would involve a detailed understanding of the geometry of the 
collision, and its consequent uncertainties. In view of these, we feel that no firm 
statements can really be made within the context of this approach, to peripheral 
collisions in heavy ion collisions. 

In conclusion, in this paper, we have used the technology of NRQCD to esti- 
mate the contribution of the rj c meson to dilepton production in the IMR region of 
heavy ion collisions, and found that it is able to satisfactorily explain the deficit, 
within the limits of the assumptions used in this study. 
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